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Abstract We study plasma flows along selected coronal loops in NOAA Ac- 
tive Region 10926, observed on 3 December 2006 with Hinode's EUV Imaging 
Spectrograph (EIS). From the shape of the loops traced on intensity images 
and the Doppler shifts measured along their length we compute their three- 
dimensional (3D) shape and plasma flow velocity using a simple geometrical 
model. This calculation was performed for loops visible in the Fe VIII 185 A, 
Fe X 184 A, Fe xn 195 A, Fe xm 202 A, and Fe xv 284 A spectral lines. In 
most cases the flow is unidirectional from one footpoint to the other but there 
are also cases of draining motions from the top of the loops to their footpoints. 
Our results indicate that the same loop may show different flow patterns when 
observed in different spectral lines, suggesting a dynamically complex rather 
than a monolithic structure. We have also carried out magnetic extrapolations 
in the linear force- free field approximation using SOHO/MDI magnctograms, 
aiming toward a first-order identification of extrapolated magnetic field lines 
corresponding to the reconstructed loops. In all cases, the best-fit extrapolated 
lines exhibit left-handed twist (a < 0), in agreement with the dominant twist of 
the region. 

Keywords: Active Regions, Structure, Velocity Field, Magnetic fields. 



1 Research Center for Astronomy and Applied Mathematics, 
Academy of Athens, 

Soranou Efcsiou 4, 11527, Athens, Greece 
email: psyntelis@phys.uoa.gr 
email: cgontik@academyofathens.gr 

2 Section of Astro-Geophysics, Department of Physics, 
University of Ioannina, 45110 Ioannina, Greece email: 
calissan@cc.uoi.gr 

3 National Observatory of Athens, 

Lofos Nymphon, Thission 11810, Athens Greece 

4 Section of Astrophysics, Astronomy and Mechanics 
Department of Physics, University of Athens, 
Pancpistimiopolis 157 84, Zografos, Greece email: 
tsingan@phys.uoa.gr 



SOLA: loop_reconstruction.tex; 6 March 2013; 14:34; p. 1 



Syntelis P. et al. 






EIS Raster 1 










LIS 


Raster 2 . 



-10001 

-1000 



Figure 1. Full disk image recorded by SOHO/MDI on 3 December 2006, at 20:51 UT. EIS 
rasters 1 and 2 field of views are represented with dark frames around AR 10926. The white 
frame shows the part of the MDI magnetogram used for the magnetic field extrapolation. 



1. Introduction 

Solar telescopes provide us with two-dimensional projections of coronal struc- 
tures on the plane of the sky. The measure of the 3D geometry of coronal loops 
is important for understanding these structures. Reconstructions of a loop's true 
shape have been attempted via a variety of methods, such as stereoscopy or mod- 



els that make assumptions for the loop shape, flow etc ( Loughhead, Wang, and Blows, 1983|| 



|Berton and Sakurai, 19 85 Asch wanden et al., 1999; Nitta , van Drie l-Gcsztc lyi, and Harra-Murnion, 1999[ ).[ 
Moreover, the STEREO space mission, with its twin telescopes in two different 
positions in space, is dedicated to the solar coronal stereoscopy by combining 
simultaneous EUV images from two different lines of sight ([Aschwanden, 2009; 



Aschwanden and Wiilser, 2011). The 3D geometry of active- region loops, com- 
puted with STEREO images, was used to constrain the magnetic field topology 
computed with non- linear force- free extrapolation methods ( |De Rosa et al., 2009[ ).[ 
To interpret the active-region Doppler maps computed with the data from in- 
struments such as Hinode/ElS and to understand the plasma flows along coronal 
loops, one also needs to know the loops' 3D geometry. IDel Zannal (|2008p . using 
Hinode/ElS observations of NOAA AR 10926, studied the behavior of the line 
of sight velocity, along loops and weak emission regions of the active region, at 
different plasma temperatures. 

In this work we analyze the same active region and we study the flows along 
selected loops for which we were able to reconstruct their 3D geometry. Our 
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Figure 2. Left: EIS intensity image of the active region in Fe XII 195 A for raster 2. Right: 
The corresponding EIS dopplcrgram. For the dopplergram, white/black colors corresponds to 
red/blue shifts. Doppler shifts are in the range of -20 to 20 km s — 1 . 



method, introduced by |Alissandrakis, Gontikakis, and Para 120081 (Paper I) also 
uses the Doppler shifts measured along coronal loops in this reconstruction. 
Based on the analysis of Paper I we extend the study to more loops and a wider 
range of formation temperatures in this work. 



2. Observations 

NOAA AR 10926 was observed on 3 December 2006, at (395", -198") from 
disk center (see Figure[l}. During that day, Hinode's EUV Imaging Spectrograph 
(EIS) performed two rasters of the AR, recording spectral lines of Fe viii 185 A, 
(10 5 - 8 K), Fe x 184 A, (10 6 K), Fe xn 195 A, (10 6 - 1 K), Fe xm 202 A, (10 6 2 K), 
and Fe XV 284 A, (10 6 3 K). These formation temperatures have been published 
by [Young et~aT\\2U07\ Raster 1 was obtained from 15:32:19 UT to 17:46:31 UT 



with a 256" x 256" field of view (FOV), while raster 2 was obtained from 
19:15:12 UT to 23:44:09 UT and covered a 512" x 256" FOV Both rasters were 
scanned from West to East (see Figure [T]) with a 30 s exposure time for each 
slit position, and a 1" spatial pixel size. We also used a timeseries of TRACE 
171 A filtergrams, recorded from 19:10:49 UT to 23:49:06 UT, with 0.5" spatial 
pixel size and a time cadence of 1 minute. Magnetic field extrapolations were 
performed using a SOHO/MDI magnetogram. The MDI full disk line-of-sight 
magnetogram (Figure [l} has a rough pixel size of 1.9" but covered the entire 
active-region ficld-of-view and was recorded on 3 Dec. 2006 at 20:51:01 UT, 
during raster 2. SOT provided a vector magnetogram with a pixel size of 0.15" 
in a FOV of 316" x 157", pointing at (447",-171"), that was taken on 4 Dec 
2006. from 10:07:00 UT to 13:03:36 UT. 



3. Data processing and loop reconstruction 

The EIS data were treated following standard procedures. After correcting for 
electron spikes, pedestal and dark current, hot pixels, slit tilt and thermal drift, 
we performed Gaussian fits on all individual line spectral profiles to compute 
Dopplcrgrams and intensity images (see Figured]). Fe viii 185 A is blended with 
Ni xvi 185.23 A (T = 10 6 4 K, |Young et at] HMD, thus a double Gaussian 
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fit was performed to estimate the effect of the blending on the dopplergrams. 
As the loops selected in Fe xn 195 A were outside the active region core, the 
Ni xvi 185.23 A line intensity is negligible, and does not influence the measured 
Dopplcr shift. 

To calibrate the zero level of the measured Doppler shifts, we averaged the 
velocities in small quiet areas near the active region observed during the second 
raster. We compared our dopplergrams with the ones published in IDel Zannal 
(2008) and found that they were in good agreement. 

For every spectral line, pointing corrections relative to He II 256 A were taken 
into account to co-align the EIS data. Furthermore, in order for the different data 
sets to be comparable, co-alignment was performed among the EIS, TRACE, 
SOT and MDI data, within the uncertainty of the MDI spatial resolution. 

From the aforementioned EIS data we manually selected, at each wavelength, 
the loops that could be distinguished from the background, along their entire 
length (Figure [3]). For each loop, we recorded their positions (x s ,y s ) on the 
intensity image as well as the Dopplcr shifts Vs z , at the same positions, from 
the Dopplergrams. We applied a boxcar average with a width of 5 to 7 at these 
line of sight velocities to reduce noise. 

The measured positions (x s ,y s ) and line-of-sight velocities Vs z along the 
loops' length, were then used as input for the computation of the loops' 3D 
structure following the methodology of Paper I. Briefly, the model assumes that 
a given loop is (i) stationary, (ii) lying in a single plane, and (hi) exhibits plasma 
flows along its length. The sought-after unknown parameter is the inclination f3 
of the plane with respect to the local solar vertical. There are two extremes for 
fa. fa, that is the line-of-sight inclination with respect to the local vertical, and 
a maximum inclination fa = fa — 90° such that the loop plane is not submerged 
below the solar surface. 

To find the unknown inclination j3 we first transform the loop's positions and 
line-of-sight velocities from the sky coordinate system (x s ,y s ,z s ) to the local 
orthogonal system (xi,yi,zt), where the (xi, ?/i)-plane is tangent to the solar 
surface and xi runs parallel to the equator, and then to the loop's system (x, z). 
Here x is defined by the loop's footpoints, z is lying on the loop's plane, and 
V(s, fa) is the velocity on the plane, along the loop's length (see Figure 3 of 
Paper I for the geometrical setup). 

We apply the method by changing the inclination within the interval [fa , fa] 
in steps of 0.1°. At each step, we monitor the quantity 

M(P) = maxAV(s,l3), where AV(s,fa = \V(s i+u /3) - V(s u fa\ is the abso- 
lute velocity difference between consecutive pixel positions, i, i+1, along the loop. 
The optimal /3 is the one minimizing M(fa. We further define an uncertainty 5 ft 
such that M{fa is smaller than, twice the minimum M(/3)-value. In effect this is 
a 2 — (t uncertainty level; see Figure HJ: for an example. Most fa values within the 
inspected range are immediately discarded either because large discontinuities 
appear in V(s,(3) along the loop or because V(s,fa obtains unrealistically large 
values. Besides the primary minimum M{fi) for j3 G [fa, fa] we sometimes obtain 
a secondary minimum for /3 ~ fa. This is also typically discarded as the loop- 
length solution in this case becomes unrealistically large, larger than the solar 
radius. 
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Figure 3. Intensity images, in various spectral lines, from both rasters, on which we indicate 
the shape of the selected loops. Loop 1 appears in all lines and in all rasters. Loop 2 can be 
seen in panels a and e, Loop 3 in panels b and f, Loop 4 in panels c and g, loops 5 and 6 in 
panels d and h. 



4. Results 

In the studied active region, during both rasters and in all recorded spectral 
lines, we were able to identify six loops having a sufficient contrast with respect 
to the background along their entire length. All of them were located in the 
south part of the active region (solar-Y coordinates less than -200"). For an 
additional four loops the method did not work successfully, failing to converge 
to a continuous velocity solution. 

Figure El shows all studied loops as they appear in intensity images of different 
spectral lines during the two rasters. Loop 1 was the most prominent, as it 
appears in both rasters and in all wavelengths shown in Figure [3J Loop 2 was 
observed during raster 1 in Fe xiii 202 A and Fe xn 195 A (Figure [3] panels a,e). 
Loop 3 was observed during raster 2 in Fe xv 284 A and Fe XIII 202 A (Figure [3] 
panels b,f). Loop 4 was observed in Fe x 184 A during both rasters (Figure [3] 
panels c,g). Finally, loops 5 and 6 were observed in Fe VIII 185 A during raster 1 
and in Fe xn 195 A during raster 2 respectively (Figure [3J panels d,h). 

4.1. Results for Loop 1 

4- 1.1. Loop reconstruction 

Figure H] shows the application of the reconstruction method for Loop 1, from 
the Fe xn 195 A image during raster 2 (see Figure [3J panel h). In the first 
row, the two first panels show the plane of the sky projection of the loop (panel 
a) and the Doppler velocity along its length (panel b), where positive values 
correspond to redshifts. The error bars in panel b are derived from the statistical 
uncertainties of the Gaussian fit. The data presented in panels (a) and (b) are the 
inputs of the reconstruction method. Then, for a wide range of inclinations /3 we 
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Figure 4. Results for Loop 1 at Fe XII 195 A. First row, panels a and b: Loop position in 
the sky plane and line-of-sight velocity as a function of the loop's length, panel c, shows M(0) 
(solid line) and z ma x{0), maximum loop altitude, (dashed line) as a function of 0. The two 
vertical lines show the range of optimum solutions while the vertical dashed line indicates the 
selected 0. Each row of four panels shows calculations for different loop inclinations, 0. These 
include: The reconstructed loop in its own plane (x, z) (panels d, h,l), the flow velocity V along 
the loop (panels e,i,m), the horizontal component V x (panels f, j, n) and vertical component 
V z (panels g, k, o) as a function of the loop length. The third row is for the optimal value of 
= —63°. The diamond indicates the top of the loop and the box indicates the west footpoint. 
The arrows in panel h) show the direction of the plasma flow. Note the discontinuities of the 
flow velocity in the second and the fourth rows (0 = —49° and = —69°). Positive line of 
sight velocities correspond to rcdshifts. 



computed the loop shape and the velocity along the loop, in the loop reference 
system (x, z). Panel c) shows the maximum velocity difference (solid line) M(/3) 
and the maximum loop altitude z max (dashed line) as a function of /3. Panel 
(c) presents a 'forest' of large M((3) values, of the order of 10 3 to 10 6 km s _1 
interrupted by a region with values of 10 to 200 km s _1 . For large j3 values, 
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M((3) has another small value region, which corresponds to unphysically large 
values for z max (~ 10 000" ). 

Each of the next three rows, has four panels, showing, from left to right, the 
reconstructed loop on its own plane, the velocity V of the flow along the loop, 
the horizontal V x and the vertical V z velocity components as a function of the 
loop length. The second and fourth row show the results for (3 = —49° and —69° 
respectively. There, the transformation from the sky plane to the loop plane 
gives a discontinuity in the computed velocity (panels e, f and g) for (3 = —49°, 
and a strong peak in velocities that does not seem physical (panels m, n and o for 
f3 = —69°), thus these results are rejected as unphysical. For the third row, where 
(3 = —63°, the velocity is continuous (panels i, j and k), and therefore the results 
are considered acceptable. The examined values of j3 range from /3a = —75° 
where the loop plane becomes tangent to the solar surface, to j3i = 14°, where 
the line-of-sight is inside the loop plane. Negative values of (3 means that the 
loop is inclined towards the south part of the active region. Let us come back 
to panel c. At /3 = —49° a square shows the M(f3) value of 4.8 x 10 4 km s _1 
which corresponds to the subtraction between the most negative values of v seen 
in panels (e), (f) and (g) at s = 19 arcsec, where V, V x and V z resemble with 
delta functions. Panel c), at /3 = —69° a square shows the value of ~ 105 km s _1 
which corresponds to the maximum difference values at s ~ 155 arcsec in panels 
(m), (n) and (o). 

The two horizontal lines show the minimum and two times the minimum of 
M(f3). These yield a range for (3 (shown between the two continuous vertical 
lines), defining an uncertainty of ±4° for the loop's inclination (3. The mean 
value of (3 in this interval is pointed by the vertical dashed line. 

4-1-2. Interpretation of flows in Loop 1 

The flow velocity V, in Figure [4ji, for (3 = —63°, is negative, with the exception 
of the west footpoint of the plot. This means that the velocity flow is towards 
the east footpoint starting from the upper part of the west leg, through the loop 
apex as indicated by the two arrows in panel h. There is however also a portion 
of the loop, near the west footpoint where we observe a downflow indicated by 
the right arrow in panel h. In panel j, V x , is negative, indicating also a motion 
towards the east footpoint. The west leg is almost vertical which is why V x ~ 0. 
In panel k, V z also shows a flow towards the east footpoint with negative/down- 
fiows, east of the loop top (diamond), and positive/up- flows west of the loop top. 
Near the west footpoint, V z changes sign again, at the same position where also 
the Doppler shift changes sign (panel b). This would mean that the plasma flows 
towards the west footpoint in this small section. However, the small Doppler shift 
values increase the relative error and makes interpretation ambiguous, for this 
small loop section. 

Figure [5] shows the results for Loop 1, for all spectral lines during raster 2. 
Each row in Figure [5] corresponds to a different spectral line and includes three 
panels. The left panel shows the loop shape on its own plane, the middle panel 
shows the measured Doppler shift along the loop and the third panel shows the 
magnitude of the computed velocity. For the middle panels, positive values mean 
redshifts and error bars arc due to Gaussian statistical uncertainties. 
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Figure 5. Results for Loop 1 for all spectral lines during Raster 2. First column: The recon- 
structed loop in its own plane (x,z). Second column: Line-of-sight velocity as a function of 
the loop's length. Third column: Velocity of the flow as a function of the reconstructed loop's 
length. 



From the first column of panels we can see that Loop 1 implies similar inclina- 
tions in all spectral lines, in the range of —63° to —60°. For the Fe vin spectral 
line, we measure redshifts of 5 to 20 km s" 1 along the loop (panel b) which for 
the given loop inclination correspond to unidirectional flow from the east toward 
the west footpoint (panel c). The results for Fe x spectral line, show that both 
the Dopplcr shift (panel e) and the computed velocity (panel f) change sign near 
the loop top; this indicates a draining motion from the loop top (indicated with 
a diamond symbol) towards the footpoints. This draining shows velocities of 
~ 40 km s -1 at both footpoints. The last two spectral lines, Fe xin 202 A and 
Fe XV 284 A show a similar dynamic behavior. The measured Doppler shifts 
are towards the blue (panels k, n) and the computed velocity (panels 1, o) shows 
unidirectional flows from west to east footpoints. This similarity indicates that 
we observe the same plasma structure in these two spectral lines. The flow of 



SOLA: loop_reconstruction.tex; 6 March 2013; 14:34; p. 8 



On the shape of active region coronal loops observed by Hinode/EIS. 



Table 1. A summary of the 3D reconstruction results for all selected loops. The Table includes: 
loop number (column 1), spectral line (column 2), raster number (column 3), j3 (column 4), 
the type of flow, which can be draining, or flow from East to West (E. to W.) or West to 
East (column 5), apparent density scale length (column 6) maximum altitude (column 7) and 
length (column 8). The last three columns, show the velocity mean values along 10% of the 
loops length starting respectively from east and the west footpoint and the mean velocity along 
10% of the loop top. The corresponding standard deviation is shown in parenthesis. 
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(°) 




(") 


(") 


(") 


East 


West 


Top 


1 Fe vin 185 A 
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-1(3) 


1 Fe xii 195 A 
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76(34) 


26(7) 


1 Fe x 184 A 


2 


-60 ±4 
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-46(6) 
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-8(3) 
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2 


-63 ±4 
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21 
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26 
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-21(3) 
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4 
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-16(1) 


-11(3) 


-16(1) 


3 Fe xv 284 A 
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-61 ±6 


W. to E. 
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55 


-11(1) 


-17(5) 


-7 (1) 


4 Fe X 184 A 


1 


-69 ± 7 


Draining 


66 


7 


73 


-33(6) 


31(4) 


-9(5) 


4 Fe X 184 A 


2 


-65 ± 7 


W. to E. 
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9 


65 


-39(4) 


-21(2) 


-6(2) 


5 Fe vm 185 A 


1 


-67 ±0.3 


E. to W. 


67 
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185 


30(25) 
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31(5) 


6 Fe xii 195 A 
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-68 ± 5 


W. to E. 


56 


24 
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-14(2) 


3(4) 


-12(2) 



Fe Xii 195 A in panel (i), has also the characteristics of a unidirectional flow 
as in panels (1) and (o) but the positive flow at the west footpoint makes this 
interpretation problematic. It is worth noting that the flow in the hot lines is in 
the opposite direction to the flow in the cooler Fe vm 185 A. Therefore, this may 
imply that Loop 1 is composed of different strands or background structures. 

4.2. Results for all loops 

In Table Q] we present a summary of the general characteristics of the flow for all 
studied loops in all the recorded spectral lines for both rasters. As seen in column 
4, all the studied loops, in the different spectral lines have a high inclination in 
the range -60° to -75° with an average inclination of -64.5° and an uncertainty 
of 3.5°. The inclination measured for Loop 1 has an average of -63.3° with an 
uncertainty of 1.8°. Loop 1 inclinations in the first raster have a small scatter, 
in agreement with the smaller error bars. The larger error on (3 is measured for 
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loop 3 (line 12) and is of ±15°. In column 5 more than half of the loop (53%) 
show unidirectional flow from West to East, 30% show draining flow from the 
loop top to the footpoints while the rest of them, all in Fe vin 185 A line, have 
flows from East to West. Three of the four cases observed in the Fe x 184 A line 
presented draining motions. 

We note that the dynamics of Loop 1 does not change qualitatively from the 
first to the second raster, except for the case of the Fe xn 195 A line. Loops 
4, 5 and 6 also change dynamics from the first to the second raster. For the 
cases with qualitatively similar loop dynamics in both rasters we can argue for a 
quasi-steady flow along them between 16:05 UT, the time when the first raster 
started scanning the loop, and 21:31 UT, the time the second raster finished, 
which is five and a half hours. Moreover, the ratios of the loop lengths (Table [IJ 
column 8) over the corresponding velocities at the footpoints (columns 9 and 
10), give timescales, 90% of which, are within 25 minutes to 2 hours. These 
values are close to the duration of the loops rastering (about 1 hour). For the 
other loops we cannot argue whether or not their dynamical behavior changes 
between the two rasters. However, as long as the loop magnetic structure is not 
modified during the raster, the applied method should give valuable results for 
the inclination /J. 

The last three columns of Table [T] show the velocity average values along 
10% of the loop's length starting, respectively, from the east footpoint the west 
footpoint and at the top of the loop, with their standard deviation in parenthesis. 
In most cases, the east footpoints have larger velocities by, factors from 1.2 to 
3. We suggest that the loops appearing in Fe xin 202 A and Fe xv 284 A 
should correspond to the same loop because the calculated dynamics is very 
similar. Columns 7 and 8 show the maximum altitude from the solar surface and 
the true loop lengths respectively. Loop 1 is apparently composed of separate 
strands that appear in different spectral lines as indicated by the different kind of 
flows found. For each spectral line, the reconstructing method gave values of the 
loop height (column 7) which vary from 23" to 26" while the loop reconstructed 
length varies from 170" to 190". These variations indicate the uncertainties of the 
method but could also show the difference between various strands composing 
the loop. 

4-2.1. Calculation of length scales and Mach numbers 

The hydrostatic pressure scale height, defined in centimeters as (A„ = 4.7 x 
10 9 (T/1MK )), ( |Aschwanden, 2005D for the spectral lines Fe VIII , 185 A, (10 5 - 8 K),| 
Fe x, 184 A, (lO 5 K), Fe xn, 195 A, (10 61 K), Fe xm, 202 A, (10 62 K), and 
Fe XV, 284 A (10 6 3 K), is 41.5", 65.7", 82.7", 104.2", and 131.2" respectively, 
and was computed using the formation temperature of each spectral line given in 
the parentheses. This scale height is still meaningful for subsonic Sows which is 
the case of our loops as it is shown in the next paragraph. Note that the maximum 
altitudes (column 7) of all measured loops are smaller than the corresponding 
pressure scale height. Thus this explains that the tops of loops are dense enough 
to be bright and detectable. 

Column 6 shows the apparent density scale, A° bs (see lAsc hwandcnl 120051 p. 
84). To calculate A° bs along the loops, we first corrected the loop intensity from 
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the background contamination. To estimate the background, we represented each 
loop in images with curvilinear grids, (jAschwanden et al.[ [T999I their Figure 6). 
The x-axis of these images represents the length along the loop so that the 
loop looks stretched and horizontal. For each point (xi oop ,yi oop ) along the loop 
in the curvilinear grid, we selected pixels lying along the line perpendicular 
to the loop {xi oop ,y) with intensities Iback{xioo P ,yback) < fIioo P {xioop,yioo P ) 
where factor / is in the range from 0.5 to 0.8. We subtract the average of 
the background pixels intensities from the corresponding loop pixel intensities 
to get the corrected intensity along the loop. We then calculate the electron 
density n e (s) along the loop as follows. For each loop we use the contribution 
functions G{T) of the corresponding spectral line from the CHIANTI software 
( |Dere et al., 1997] ) . The contribution function for all spectral lines was computed 
assuming a IMazzotta et al.\ 119981 ionization fraction and a hybrid abundance 
for iron ( |Fludra and Schmelz, 1999) ). We assume that the loop width along the 
line of sight is constant and of w — 2". The electron density is calculated as 
n e (s) = \J w q^P — -■ T max is the spectral lines formation temperature, maximiz- 
ing the contribution function and s is the length along the loop. We performed a 
fit with an exponential function to the electron density function n e (s) computed 
along the loop. The density scale length A° bs is the derived length from the 
exponential fit. The resulting relative difference between A° &s and X p /cos(f3) is 
in the range of 0.01 to 0.8 with an average of 0.46. 

The sound velocity, expressed as v s — 0.151 x VT km s _1 , was computed 
for the different spectral lines formation temperatures to be in the range of 
120 km s _1 to 213 km s _1 . Therefore, the corresponding Mach numbers, are 
smaller than one with the exception of the west footpoint of loop 5 (row 15 
in Table [1]) which is supersonic. For the Fe vm , 185 A, Fe x , 184 A, Fe xn , 
195 A, Fe xiii , 202 A and the Fe xv , 284 A lines, the mean Mach numbers 
computed from the corresponding velocities in columns 9 and 10 (Table [1} are 
0.6, 0.3, 0.2, 0.2 and 0.1 respectively and shows that, the larger the line formation 
temperature the smaller the Mach number. 

4-2.2. Influence of background correction on Doppler shifts 

In order to check the influence of the background on our results, we also per- 
formed the loop reconstruction, using loop Doppler shifts calculated from spec- 
tral profiles from which we subtracted a background spectral profile ( |Gontikakis et al., 2005[ ) 
We got better results when the background points were selected along the same 
slit position as the corresponding loop points. We obtained similar inclinations 
within ~ 3° compared to the uncorrected case. Qualitatively the flow (draining 
or unidirectional flow) is the same for nine cases (rows 1, 2, 3, 6, 7, 10, 14, 15, 16 
of Table [l}. The calculated flow, in most of the other cases diverged at the west 
leg of the loops. The west leg of the loops is almost parallel to the slit so that 
our method to select appropriate background points fails. 
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5. Comparing reconstructed with extrapolated coronal-loop shapes 

As an independent test of our 3D reconstructions we used linear force-free 
magnetic field extrapolations as described by lAlissandrakisI (|198ip to deter- 
mine whether extrapolated field lines (viewed as loops) that model the recon- 
structed loops can be found. That active-region magnetic fields are far from a 
linear force-free configuration (and, in fact, far from a force-free field config- 
uration in general, at least at photospheric heights) is well known (see, e.g., 
IMetcalf et~al[ [T9951 |Georgoulis and LaBohte] [2004)1 . However, we arc seeking 



single a- values to fit single loops, with different values corresponding to different 
loops. In this approach we have better control over the fitting process, that 
cannot be achieved by a fully-fledged model-dependent nonlinear force-free field 
extrapolation (see, e.g., |Schrijver et ~aT\ l2006UMetcalf et aU [20 08). To compare 
the observed and the extrapolated loops we employ the distance modulus C;, 
introduced by ( |Wiegelmann and Neukirch, 2002| ), namely, 

Cf = —J / \Hloop — R-line\dl , (1) 

L loop J 

where R/ D op, R-Hne are respective vector positions along the reconstructed loop 
and extrapolated field line, respectively, and Li oop is the loop's length. In- 
tegration occurs along the reconstructed loop's length with respective points 
considered as points having the same location (in terms of normalized distance 
from footpoints) along the loop. The number of points taken vary with the loops' 
length: for the smallest loop (loop 3) we used 51 points while for the largest (loops 
1 and 5) we used up to 215 points. We only compared loops and extrapolated 
field lines with footpoints lying within a square with a 30"-linear size. 

As lower boundaries for the extrapolations we use a Level 1.8 magnetogram 
from the full-disk Michelson Doppler Imager (MDI; IScherrer et al.\ [T995) on- 
board SOHO. The course of action is to, first, crop the MDI magnetogram to 
select the desired photospheric patch and then use the cosine-corrected line-of- 
sight magnetic field component on the image (sky) plane. We choose this action 
because de-projecting for the local, heliographic plane would introduce unfore- 
seen uncertainties when similarly de-projecting the EIS images. Point taken, we 
understand that each approach has weaknesses - this is an issue that we plan to 
investigate thoroughly in future studies. 

The linear size of the selected patch is L — 183 Mm. This determines the ex- 
tremes of the |a|-range to be investigated fe.g..lAlissandr akislll981l [Green et al\ 
I2002p . by the formula |a maa; | = (2n/L). Hence, we scan an a-range between 
-0.0343 Mm" 1 and 0.0343 Mm -1 , with a step-size equal to 0.006 Mm" 1 . 

An example loop fitting, corresponding to Loop 1 at Fc xn 195 A is shown 
in Figure El The minimum C;-value achieved in this case (Figure [B^; also shown 
in Table 2) is C/ ~ 0.2. Figure [5] also includes the best-fit extrapolated field line 
(thin solid curve) and a potential-field line (dotted curve) for reference. Notice 
that, despite the proximity between the best-fit line and the reconstructed loop, 
the match is not perfect. In particular, the best-fit line shows smaller inclination 
than the loop and larger distances between the footpoint locations as compared 
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Figure 6. Fitting Loop 1 at Fe XII 195 A during raster 1: (a) the C; modulus along the loop 
as a function of a. A minimum C; is achieved for a ~ —0.009 Mm~ . (b, c) Two different 
views of the inferred loop (thick solid curve) , accompanied by a potential- field extrapolated line 
(a = 0; dotted curve) and the "best-fit" extrapolated field line (a = —0.009 Mm -1 ; thin solid 
curve), (d) The same system, seen from above and projected on the respective SOHO/MDI 
magnetogram. 



to the distances between the two (loop and field-line) apices. In addition, the 
C/-curve is quite shallow (as also reported by |Wiegelmann and Neukirch I2002[) 
so different a-values give rise to very similar extrapolated lines. In situations 
like this it is challenging to assign an uncertainty to each best-fit value. For the 
purpose of providing an error margin in the best-fit a-values, however, we have 
selected a 10%-margin in C; enclosing the minimum and calculated the standard 
deviations of a- values within this margin. For the example of Figure [5] we find 
a = —0.009 ± 0.006 Mm -1 . Even given the uncertainty, therefore, Loop 1 at 
Fe XII 195 A is better fitted by a negative-a field line. 

From the 17 loop cases listed in Table 1 the fitting method gave a minimum 
(albeit shallow) C\ in 15 cases. The results are summarized in Table 2. From 
them, we notice that, first, all best-fit a-values, including uncertainties, are neg- 
ative, indicating left-handed twist in the best-fit extrapolated field lines. Second, 
the distance Z between the solar surface and the reconstructed loop apices is wcll- 



S0LA: loop_reconstruction.tex; 6 March 2013; 14:34; p. 13 



Syntelis P. et al. 



Table 2. Summary of the results obtained from the magnetic field extrapolation in comparison 
to the results of the 3D reconstruction. Columns 4 and 5 are the reconstruction results as in 
Table ^ Columns 6 to 8 are the extrapolation results for the best fitting magnetic line and 
column 9 is the calculated best a with its error bar. 



loop 


spectral 


rast. 


P 


Z 


P 


Z 


Ci 


a ± 5a 


num. 


line 


num. 


(°) 


(") 


mag.(°) 


mag.(") 




(Mm- 1 ) 


1 


— 

Fc vin 185 A 


1 


-64 


26 


-55 


26 


0.23 


-0.013±0.006 




Fe x 184 A 


1 


-65 


25 


-52 


25 


0.23 


-0.012±0.007 


) 


Fe xii 195 A 


1 


-65 


23 


-51 


24 


0.20 


-0.009±0.006 




Fe xiii 202 A 


1 


-65 


23 


-51 


24 


0.24 


-0.009±0.007 




Fc vin 185 A 


2 


-61 


26 


-50 


25 


0.13 


-0.006±0.005 




Fe x 184 A 


2 


-60 


25 


-50 


22 


0.15 


-0.008±0.004 




Fe xii 195 A 


2 


-63 


24 


-50 


25 


0.15 


-0.006±0.004 




Fe xiii 202 A 


2 


-63 


25 


-50 


25 


0.15 


-0.006±0.004 




Fe xv 284 A 


2 


-64 


26 


-48 


30 


0.14 


-0.003±0.004 


2 


Fe xii 195 A 


1 


-64 


15 


-44 


13 


0.23 


-0.010±0.009 


2 


Fe XIII 202 A 


1 


-64 


15 


-47 


15 


0.23 


-0.015±0.009 


3 


Fe XIII 202 A 


2 


-75 


1 


-49 


11 


0.19 


-0.010±0.004 


3 


Fe xv 284 A 


2 


-61 


6 


-50 


8 


0.22 


-0.004±0.006 


5 


Fe vin 185 A 


1 


-67 


27 


-51 


27 


0.24 


-0.007±0.007 


6 


Fe xii 195 A 


2 


-68 


7 


-54 


12 


0.12 


-0.007±0.002 



reproduced by the fitting. Excluding Loop 3 at Fe xiii 202 A (7" of difference) 
and Loop 6 at Fe xn 195 A (5" of differ ence), the mean amplitude difference 
in Z between loops and extrapolated field lines is (1.2 ± 1.3)", corresponding to 
~ 5.5% of the mean loop length. When it comes to loop inclinations, however, 
extrapolation gives consistently lower inclinations. The mean inclination differ- 
ence between all reconstructed loops and extrapolated field lines in Table 2 is 
14.5° ± 4.3°, or (22.5 ± 5)% of the reconstructed loop length. The reason(s) for 
this systematic difference need to be investigated further. Generally speaking, 
however, more than one effect may be contributing to this: besides the limited 
validity of the force-free approximation in the photospheric magnetic fields used 
as extrapolation boundaries that, however, is not clear whether should yield a 
systematic lower inclination, (i) the central assumption of loops lying in a single 
plane may overestimate the loops' inclination as the most highly projected parts 
of each loop weigh heavier in the reconstruction, and (ii) the extrapolations using 
the image-plane magnetic field include some projections that may well give rise 
to an underestimation of the extrapolated field lines' inclination. 

Summarizing, it is indeed a nontrivial exercise to compare reconstructed 
loops with extrapolated field lines. Previous works (Wicgclm ann and Neukirch] 
2002] IDemoulin et alX 120021 [Carcedo et al. . 2003] and others) also reach similar 
conclusions. Why is this first-order comparison attempted here, then? A rather 
straightforward answer is because, despite difficulties, (a) we closely match the 
loop apices' height from the solar surface, and (b) all best-fit lines correspond 
to negative a-values, in statistical agreement with the best-fit a-value obtained 
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for the entire NOAA AR 10926 (a = -0.049 ± 0.006 Mm' 1 ), as calculated by 
the method described in |Georgoulis and LaBonte| 120071 using the SOT mag- 
netogram. This implies that the linear force-free method applied separately to 
each loop consistently yields a dominant left-hand twist in the AR, in agree- 
ment with the cruder linear force-free approximation for the entire AR. Even 
the systematic lower inclination for extrapolated field lines is an effect that is 
interesting to investigate further, as one may be able to assess the validity of 
the loop reconstruction method vis d vis various effects in the extrapolation. We 
intend to continue and extend these comparisons to clarify these issues. 

6. Summary and Conclusions 

We reconstructed the 3D geometry of six loops observed by the Hinode/FilS 
spectrograph in five spectral lines during two rasters of the instrument, a total 
of 17 cases. All loops correspond to NOAA AR 10926. All reconstructed loops 
have large inclinations with respect to the vertical, in the range of —60° to —75° 
and they are all projected to the South of the bright central part of the active 
region. Moreover, due to pressure scale height effects, large inclinations lead to 
denser and brighter plasma near the loop top, hence the entire loop stands clearly 
above the background level. The flows calculated for all loops were subsonic in all 
spectral lines that validates the loops description using their hydrostatic scale. 

The best-studied loop is Loop 1, successfully reconstructed in nine images. 
We find that this loop is not a monolithic structure, as the flows deduced in 
different spectral lines vary from unidirectonal flow from East to West in the 
low-temperature Fe VIII line, to draining motion from the top to the footpoints 
in the intermediate-temperature Fe x line, to unidirectional flow from West to 
East in the high-temperature lines. This being said, the computed inclination is 
very similar in all images, namely 63° ± 3°, that strengthens our assessment of 
internal structure in the loop in different temperature ranges. 

An independent comparison between the reconstructed loops and extrapo- 
lated field lines by means of a linear force-free extrapolation implemented on 
a case-by-case basis give results that call for additional investigation. We have 
been able to closely model the height of the loops' apices and we have noticed 
best-fit lines with a consistent (left-handed in this case) twist in the AR, but 
extrapolated field lines are consistently less inclined than reconstructed loops. 
This discrepancy may be due to drawbacks in the reconstruction method, weak- 
nesses in the extrapolation and/or the extrapolated boundary, or a combination 
of both. Aiming to validate our technique as reliably as possible, we intend 
to carry out similar investigations in future works, relying on larger statistical 
samples of reconstructed coronal loops. 
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